Bouncing of charged droplets: An explanation using mean curvature flow 
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Two oppositely charged droplets of (say) water in e.g. oil or air will tend to drift together 
under the influence of their charges. As they make contact, one might expect them to coalesce and 
form one large droplet, and this indeed happens when the charge difference is sufficiently small. 
However, Ristenpart et al discovered a remarkable physical phenomenon whereby for large enough 
charge differentials, the droplets bounce off each other as they make contact. Explanations based 
on minimisation of area under a volume constraint have been proposed based on the premise that 
consideration of surface energy cannot be sufficient. However, in this note we explain that on the 
contrary, the bouncing phenomenon can be completely explained in terms of energy, including an 
accurate prediction of the threshold charge differential between coalescence and bouncing. 
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Droplet motion induced by electrical charges has been 
extensively studied at least since the 19th century (see 
e.g. d]; for an overview and list of references, see [2]). 
Such a process occurs in a wide variety of applications 
such as storm cloud formation, commercial ink-jet print- 
ing, petroleum and vegetable oil dehydration, electro- 
spray ionization for use in mass spectrometry, electrowet- 
ting and lab-on-a-chip manipulations (see e.g. [5Ull)). 
most of which are only partially understood. 

In this note, we are interested in the physics of two 
droplets of a conductive fluid such as water, that move 
within a somewhat electrically insulating immiscible fluid 
such as oil. The precise sizes of the droplets are not im- 
portant, and we will initially consider a column of water 
as one of the droplets, supporting a layer of oil above, 
and a charged droplet of water within the oil (see the 
left-hand side of Figure [T]). Under the influence of the 
charge difference the water droplet will then approach 
the water column until they touch, at which point one 
might expect the droplet to be absorbed into the water 
column. While this absorption does indeed occur when 
the droplet 's charge is small enough, it is a remarkable 
discovery of Ristenpart et al pjj that highly charged 
droplets tend to bounce off the water column. This is 
illustrated in the experimental observations of Risten- 
part et al that can be found in Figure [T] in which the 
charge arises from a voltage gradient introduced across 
the oil. Various explanations have been given for this 
phenomenon 2^, starting with the premise that dy- 
namics driven by surface tension would always lead to 
coalescence. However, the purpose of this note is to ex- 
plain that on the contrary, the bouncing/coalescing can 
be explained purely by postulating that after touching, 
when the charge difference disappears, the droplets move 
solely to reduce their surface area/energy as quickly as 
possible. We are then able to invoke some deep ideas of 
the so-called mean curvature flow as we explain below. 

Our starting point is the well-understood principle that 
as the droplets approach each other, they deform into a 



conical shape (cf. Taylor cones under the in- 

fluence of the differing charges, and as they touch, the 
droplets will look locally like a double cone with a thin 
connecting bridge, as illustrated in Figure [2] Experiment 
shows [12] - see Figure |4] ~ that the cone angle 9 in the 
regime we are considering (see Figure [3| is proportional 
to the field strength, and that there is a critical thresh- 
old angle Oq below which one observes coalescence, and 
above which one observes bouncing, independently of the 
choice of fluids. The diagram in Figure |4] suggests that 
the critical contact angle 6c should be about 27°. 

Our goal in this note is to survey the theory of mean 
curvature flow from cones, which governs how a cone can 
evolve in order to reduce its surface area as quickly as 
possible, and show that the behaviour of this flow un- 
dergoes a bifurcation at a critical cone angle of 24°. We 
then indicate how abstract arguments prove rigorously 
that evolutions of the touching droplets must show co- 
alescence for cones angles below the critical angle, and 
must show bouncing for angles above the critical angle. 

It is a feature of our approach that we need make virtu- 
ally no physical assumptions or symmetry assumptions. 
We know from 2j that the charges are equalised at the 
moment of touching, and thus we postulate that after 
touching, the droplets evolve in order to reduce their sur- 
face area as quickly as possible. However, this evolution 
is occurring only locally in time and space, so it is irrel- 
evant that global minimisation of surface area would al- 
ways result in coalescence. Ultimately, all the subtlety of 
the distinction between coalescence and bouncing arises 
from the mathematical theory of mean curvature flow. 

Consider a one-parameter family of surfaces Mt in 
that can be viewed as the images of immersions Ft : M — > 
M.^ from a fixed surface M . When the area At of these 
surfaces is finite, it evolves by 

where the normal vector field H on the surface is called 
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Figure 1. Coalescence and bouncing of water in oil 
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A high-voltage (HV) provides the electric field of strength E, The top row of images shows coalescence (E < Ec ) whereas the 
bottom row shows bouncing (E > Ec ) of the water droplet. Reprinted by permission from Macmillan Publishers Ltd: Nature, 
WD Ristenpart et al. 461, 377-380, ©2009, doi:10.1038/nature08294. 



Figure 2. Conical structure 



Figure 3. Conical structure 2 




As the droplets touch, they are locally conical, and form a 

thin connecting bridge. Reprinted by permission from 
Macmillan Publishers Ltd: Nature, WD Ristenpart et al. 
461, 377-380, ©2009, doi:10.1038/nature0829l 



the mean curvature (see [3T1 Appendix A] for a direct 
definition). The evolution of the surface that reduces 
the area as quickly as possible (more precisely the L^- 
gradient flow) is then the much-studied mean curvature 
flow (MCF) |2T] defined by the nonlinear PDE 



dt 



Ff=H. 



(2) 



The comparison principle |21| tells us that two different 
solutions that are disjoint initially at time t = 0, will 
remain disjoint for later times t > 0. 

We locally model the touching fluid droplets as follows 
(Figure [5] shows the construction). We rotate the func- 
tion ue{xi) := \xi\/ tan 6* around the xi axis in to get 
a double cone Mg, and model the fluid bridge with an 
arbitrary one-sheeted surface Mo that smooths out Mq 
as in Figure [5) 

We can solve the mean curvature flow starting with this 
Mq, and it will either flow for all time (corresponding to 
coalescence) or will develop a singularity at finite time 




The double cone angle is proportional to the field strength 
and determines the behaviour of the system. 



as the thin neck shrinks to nothing (meaning the bridge 
between the droplets breaks apart and we get bouncing) 

To understand which of these occurs for each 0, we 
consider self-similar solutions of the mean curvature flow 
starting at the double cone of angle 6*, i.e. solutions Mt of 
the form y/iMf where Mf asymptotically approaches the 
double cone Mg . In general, there will be many solutions 
of this form with the same cone as the initial surface, 
but if we impose the ansatz that the solution should be 
symmetric under rotations about the xi-axis, and under 
reflections xi —xi, then all solutions can be classified 
up to the solution of a simple ODE PH US] . It turns out 
that for smaller than some critical value, which can 
be computed numerically to be 6% ~ 23.96°, there 
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Figure 4. Double cone angle and potential 
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The contact angle 6 of the double cone for water droplets in 
air with applied potential Aip. Open symbols denote 
coalescence, filled symbols denote bouncing. Reprinted by 
permission from Macmillan Publishers Ltd: Nature, WD 
Ristenpart et al. 461, 377-380, ©2009, 
doi:10.1038/nature08294. 



Figure 5. Double cone smoothings 




We smooth out the double cone Mq in order to obtain a 
rotationally symmetric surface Mq, locally modeling the 
touching fluid droplets. 



are at least two such solutions that are one-sheeted (see 
Figure [6]d), but as 9 increases to 6'^, these two solutions 
merge into one, a bifurcation occurs, and for 9 > 6'^, 
there is no such one-sheeted solution, only a two-sheeted 
solution which corresponds to bouncing (see Figure [6^). 

Meanwhile, a simple application of a suitable compar- 
ison principle (see |20j ) proves that for cone angles less 
than any solution of MCF starting at the smoothed 
cone Mq must evolve as a one-sheeted solution that cor- 
responds to coalescence (see Figure [6]3.). A more sub- 
tle argument [24] [25] tells us that for cone angles greater 
than 6'^, any solution of MCF starting at Mq must evolve 
as a two-sheeted solution that corresponds to bouncing 
(see Figure[6^.). For this latter case, one imagines an en- 
velope of all possible solutions emanating from the cone. 



Figure 6. One-sheeted and two-sheeted evolutions of double 
cones 




a. The two-sheeted evolution of exists for any angle 9. It 
is the unique evolution of Mo for 9 > 9q and corresponds to 
the droplets bouncing off each other. 




b. For 9 < 9q at least two one-sheeted evolutions of Mq 
exist. These can be used as pinching barriers and therefore 
imply coalescence of the droplets. 

whether rotationally symmetric or not, and argues that 
this envelope itself must be a new, regular solution that 
is now necessarily rotationally symmetric. If there exists 
any non-bouncing solution starting at a cone, then this 
envelope solution must be a one-sheeted solution of the 
type classified above, which does not exist for cone angles 
greater than e*c pill^. 

Our theory therefore predicts a transition between co- 
alescence and bouncing at a critical angle which is within 
a few degrees of the experimentally observed value. 

Finally we want to compare our approach to the one 
from 2] . There it is assumed that the bridge between the 
touching droplets minimizes area under a volume con- 
straint which leads to constant mean curvature surfaces 
of revolution (Delaunay surfaces). Similarly to our def- 
inition of double cone smoothings such a bridge is then 
fitted to a linear double cone and it is assumed that the 
sign of the mean curvature H (respectively capillary pres- 
sure) determines the behaviour. The construction is not 
differentiable and therefore there are several ways to fit 
a given Delaunay surface to a given double cone, leading 
to different predictions of 9c- The way the fitting is done 
in |2] leads to a prediction of 9c = 30.8°. 

Conclusion: We propose a new model to explain the 
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remarkable phenomenon of bouncing and coalescence of 
charged fluid droplets. Our model assumes that the sys- 
tem after touching moves according to the mean curva- 
ture flow, which means its area decreases as fast as possi- 
ble. Analyzing the flow in this setting leads to a predic- 
tion of about 24° for the critical angle, which is in good 
agreement with experiments. Therefore minimization of 
energy can, contrary to general belief, explain the phe- 
nomenon. One advantage of our approach compared to 
existing ones is that we do not make strong assumptions 
on the precise shape of the bridge between the touching 
fluid droplets. 
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